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We perform ID particle-in-cell (PIC) simulations using OSIRIS, which model a short-duration (~ bOOuj^ 1 
FWHM) scattered light seed pulse in the presence of a constant counter-propagating pump laser with an 
intensity far below the absolute instability threshold. The seed undergoes linear convective Raman ampli- 
fication and dominates over fluctuations due to particle discreteness. Our simulation results are in good 
agreement with results from a coupled-mode solver when we take into account special relativity and the use 
of finite size PIC simulation particles. We present linear gain spectra including both effects. Extending the 
PIC simulations past when the seed exits the simulation domain reveals bursts of large-amplitude scattering 
in many cases, which does not occur in simulations without the seed pulse. These bursts can have amplitudes 
several times greater than the amplified seed pulse, and we demonstrate that this large- amplitude scattering is 
the result of kinetic inflation by examining trapped particle orbits. This large-amplitude scattering is caused 
by the seed modifying the distribution function earlier in the simulation. We perform some simulations with 
longer duration seeds, which lead to parts of the seeds undergoing kinetic inflation and reaching amplitudes 
several times more than the steady-state linear theory results. Simulations with continuous seeds demonstrate 
that the onset of inflation depends on seed wavelength and incident intensity, and we observe oscillations in 
the reflectivity at a frequency equal to the difference between the seed frequency and the frequency at which 
the inflationary SRS grows. 
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I. INTRODUCTION 

Backward stimulated Raman scatterimJ^ (BSRS) in 
plasmas, in which an incident light wave in a plasma 
decays into a backward-propagating light wave and a 
forward-propagating plasma wave, has been a subject of 
much study, in large part because it scatters light away 
from the target in inertial confinement fusion (ICF)P^ 
Early research focused on relatively high intensities, 
where growth was in the weakly damped convective and 
absolutely unstable regimes, and saturation occurred due 
to wave-breaking and/or pump depletion) 5 -^ The compe- 
tition between back-, forward-, and side-scatter was also 
investigated!^ As laser and plasma parameters for ICF 
evolved, research in BSRS shifted to the strongly damped 
regime. In modern experiments, SRS typically occurs at 
densities and temperatures for which kX^e ^ 0-3, where 
Landau damping is significant. In this kinetic regime, 
the measured BSRS reflectivities can greatly exceed the 
values from linear theory calculations, e.g. in the single- 
hot-spot experiments of Ref. [9j 3D paraxial-envelope 
simulations with linear damping, however, correctly mod- 
eled the intensity threshold for SRS in experiments with 



a smoothed, multi-speckle beam in a pre-formed uniform 
hohlraum plasmaP^ A process called kinetic inflation was 
proposed to explain the single- hot-spot results JHHUl j n kj_ 
netic inflation, a small-amplitude plasma wave excited in 
the strongly damped convectively unstable regime can 
trap electrons, modifying the distribution function so 
that the kinetic damping of the plasma wave is greatly re- 
duced or vanishes ! 14 ! 15 ! Therefore, for the same incident 
laser intensity, SRS can then transition to the weakly 
damped or absolutely unstable regime. 16 There has also 
been recent work on how BSRS in the kinetic regime can 
saturate due to nonlinear frequency shift or related 
trapped-particle instabilities^ caused by electron trap- 
ping. Recent research has demonstrated the importance 
of the propagation and evolution of plasma wave pack- 
ets, including how the reflected light can occur in bursts 
spaced proportionally to the inverse of the nonlinear fre- 
quency shift P2I The latest research has also demonstrated 
that hot electrons and back- and side-scattered SRS pro- 
duced by one speckle interact with neighboring speckles, 
causing the speckles to self-organize and produce coher- 
ent bursts of SRSP^ Until recently, little work has ex- 
plored the possibility of scattered light, plasma waves, or 
the resulting changes to the electron distribution in one 
region of space or time enhancing SRS at different times 
or locations. 
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In this paper, we make a detailed comparison of the 
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linear amplification of a well defined counter-propagating 
seed pulse using coupled- mode theory and OSIRIS PKp3 
simulations. We then explore how this seed pulse can 
trigger large reflectivities after it has left the plasma us- 
ing OSIRIS simulations. We consider situations in which 
no BSRS occurs with only the pump (no seed is used). 
The seed intensity and pulse length are varied. For short 
seed pulses (~ SOOuJq 1 FWHM), the seed pulse is linearly 
amplified as it transits the box. This amplification agrees 
with linear theory when appropriately modified to take 
into account special r elativ ity and the use of finite-size 
particles in PIC codes j 21 * 22 lsuch as OSIRIS. We simulate 
ICF-relevant laser and plasma conditions, and demon- 
strate that special relativity increases the linear gain and 
shifts down the scattered light wavelength. We find, for 
these short seed pulses, that kinetic inflation occurs after 
the seed pulse leaves the box. The timing and ampli- 
tude of the first peak in reflectivity after the seed pulse 
depends on the duration and intensity of the seed pulse. 
We examine the trapped particles to verify that kinetic 
inflation is occurring and that the bounce period is con- 
sistent with the Langmuir wave (which we also call the 
plasma wave) amplitude. 



For longer seed pulses, the inflationary burst of scat- 
tered light overlaps with the seed. Under these condi- 
tions, the measured gain of the seed can reach several 
times the steady-state linear gain value when the seed 
wavelength is near the peak of the gain spectrum. We 
also examine the onset of inflationary scattering and the 
bursty nature of BSRS using continuous seeds. Non- 
resonant seeds, which are not at the peak of the linear 
gain curve, require higher incident intensity to cause in- 
flation. We also observe that, when the seed frequency 
is non-resonant, the reflected light is modulated with a 
period inversely proportional to the difference between 
the seed and resonant frequencies. (In this paper, we use 
the term "resonance," where 1 + \r — for electrostatic 
waves, and "peak gain" interchangeably.) After inflation 
sets in, the measured gain of the seed decreases with in- 
cident seed intensity due to pump depletion. 



The paper is outlined as follows. We present the sim- 
ulation geometry and plasma conditions in Section [ll[ 
and discuss the linear theory of convective BSRS gain 
and its relativistic and PIC modifications in Section Hill 
We describe in Section |IV| a subtraction technique that 
we use in our data analysis. In Section |V| we discuss 
the amplification of short-duration seed pulses in PIC 
simulations and as calculated by a coupled- mode solver. 
Section PVTl covers our observation of kinetic inflation that 
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occurs after the seed pulse passes, and Section VII 




our measurements of kinetic inflation using longer dura- 
tion seeds. Finally, we discuss the onset of inflation with 
continuous seeds in Section IVlIII and conclude in Section 

El 
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FIG. 1. (Color online) The geometry of the OSIRIS simula- 
tions. 



II. GEOMETRY AND PLASMA CONDITIONS 

Throughout this paper, we use normalized units to de- 
scribe the OSIRIS simulations. To make connection to 
parameters of interest for ICF, we assume that the inci- 
dent pump has a wavelength of Ao =351nm. When we 
note quantities in physical units, they correspond to this 
pump wavelength. Our formulas and other quantities are 
given in CGS units while the temperature is often given 
in units of eV. 

Figure [I] depicts the usual simulation geometry. The 
pump laser (cjo? &o) is incident from the left with a nor- 
malized electric field amplitude Eoe/m e cujQ = 3.68 x 10 -3 
(in physical units, Iq = 1.5 x 10 14 W/cm 2 ). The seed 
(cJi,fci) enters the simulation box from the right and 
beats with the pump, inducing an electrostatic wave 
(u)2, £2) m the plasma, which travels to the right. Light 
from the pump scatters off the plasma wave and ampli- 
fies the seed. For the simulation parameters, this process 
is a convective instability. 

Our simulation plasma has fixed ions and we do not 
add collisions. We simulate densities and temperatures 
for which SRS is expected to occur for ignition experi- 
ments at the National Ignition Facility!^ The plasma has 
a uniform density n = 0.12n c (~ 1.1 x 10 21 cm -3 ), where 
n c is the critical density for the pump laser. The electron 
thermal speed is v th = 0.0699c (T e = 2.5 keV). The box 
has length L = 1790c/u;o (100 /xm) and 8192 cells, giving 
a cell width of 1.1 A^, where A^ is the Debye length. We 
use 16,384 particles per cell and quadratic splines for the 
particle shape to reduce spurious noise due to aliasing. 
The particle boundary conditions are thermalizing and 
we use perfectly matched layers 2 -^ for the field boundary 
conditions. The particle pusher is relativist ically correct 
in all the simulations except where we state otherwise. 

The seed intensities we use in our simulations are 
far above the background electromagnetic noise in ICF, 
which for BSRS is primarily Thomson scattering! 2 ^ For 
the plasma conditions used in our simulations and a typ- 
ical NIF quad (four laser beams arranged in a square) 
with effective F-number^ 3 -! of 8 as the pump, the ef- 
fective Thomson scattering seed within the FWHM of 
the peak relativistic gain (discussed below) is 1.6 x 10 4 
W/cm 2 = 1.1 x 10 -10 io- SRS growth from such noise 
is initially linear, and enhanced over plane- wave growth 
by intense speckles in a phase-plate-smoothed beam. In 
down-stream regions, this light acts as an SRS seed far 
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above thermal noise, and may reach amplitudes where 
kinetic effects are significant. We choose seed intensities 
to induce such effects. 



III. LINEAR THEORY OF CONVECTIVE BSRS GAIN 
A. Non-Relativistic Theory 

Our PIC simulations never reach a steady state. How- 
ever, linear theory states that the reflectivity will quickly 
reach a steady state when we use a continuous scattered 
light wave seed. In the convective steady state, the seed 
intensity is amplified by a factor of e G by the time it 
exits the box. G is the linear intensity gain exponent, 
commonly called "the gain," and we present an equation 
for it in the strong damping limit. Here, we summarize 
the results from a detailed derivation for the steady-state 
gain from Ref. [25l 

Given the pump (gjo? &o) and the seed (cji, fci), we 
calculate the plasma wave (co>2, fe) using the matching 
conditions, 



and 



ko = ki + k 2 . 



(la) 



(lb) 



We need U2 and k 2 for the plasma susceptibility when we 
calculate the gain. 

Let the seed intensity be denoted by Ii(z), where z = 
(0, L) is the (left, right) edge of the box. Then, 

h(z) = e G '^h(L), (2) 
where Gi(z) is the linear intensity gain exponent, 



G t (z)= f T l {z')I {z')dz' . 

J Z 



(3) 



We neglect pump depletion and light wave damping, so 
Iq is constant, which leads to 



Ti = r s Im 



(1 + X/) 



(4) 



where the subscripts e and I denote the electron and ion 
species, respectively, Xj 1S the collisionless susceptibility 
for species j, £(^2,^2) = 1 + Y^j Xj(^2,^2) is the plasma 
dielectric function, and 
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where r e = e 2 /m e c 2 is the classical electron radius. For 
a Maxwellian velocity distribution, Xj is given by 



Xj(k 2 ,uj 2 ) 
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where CJpj is the plasma frequency of species j, VTj = 
^Tj/rrij is the thermal speed of species j, and Z'(s) = 
dZ/ds. Z'(s) must be calculated numerically, and is typ- 
ically found by first computing Z(s), the plasma disper- 
sion function.^ z'(s) = -2sZ(s) - 2, with 



Z(s) = iV^e~ s [1 + erf (is) 



(7) 



Xj — >• —{(jO P j/uo2) 2 as Tfij — >• 00, so we can set 1 + Xi ~^ 1 
everywhere (recall we use fixed ions in the OSIRIS sim- 
ulations). In particular, e — 1 + Xe and 



Ti = rjm 



(8) 



We further simplify Eq. [3] since we are dealing with a 
uniform plasma. I\ is constant, so the gain is given by 



G l (z)=T 1 I (L-z). 
We also define an amplitude gain rate, go, as 
Ti/o Xi 

go = — oc 



2 -(1 + Xp )2 + X 2' 



0) 



(10) 



where Xr and Xi are the real and imaginary parts of Xi 
respectively. We plot the theoretical gain spectrum for 
the conditions of interest in Figure [2] as dash-dotted lines. 

Equations |2][5] are valid in the strong damping limit. 
This limit applies when \v g 2da,2/dx\ <C |^2&2|, where v g 2 
is the plasma wave group velocity, CL2 is th e plasma wave 
action amplitude (defined in Section |VI3 ), and V2 is the 



Landau damping ' rate. In a homogeneous plasma in 
the convective steady- state, which is where our gain cal- 
culation applies, this condition is go <C ^2/^2- Work- 
ing at the peak of the non-relativist ic gain curve, we 
have a spatial gain rate of g% r = 3.28 x 10~ a uj /c and 
l/ 2/v g 2 = 0.061 Iujo/c. Therefore, we are in the strong 
damping limit. 



B. Relativistic Modification 

We now explore the impact of special relativity on lin- 
ear gain. Estabrook and KrueJ^ included an analysis of 
SRS for temperatures for which relativistic corrections 
are important, and performed 1.5D PIC simulations of 
laser and plasma conditions where the plasma wave is 
weakly damped. They found that non-relativistic linear 
theory does not adequately describe the wavenumber of 
the fastest growing mode in high-temperature (~64 keV) 
plasmas, but taking into account the effective (reduced) 
plasma frequency and corresponding density due to spe- 
cial relativity brings theory and simulation into better 
agreement. More recently, Bergman and Eliasson de- 
rived a fully relativistic expression for the unmagnetized 
plasma dielectric function,^ and Bers et. al. derived ap- 
proximate expressions relevant to current and near-future 
deuterium-tritium fusion plasmas P^Palastro et. al. have 
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also derived a fully relativistic description of Thomson 
scattering 

In this subsection, we simply make some heuristic 
changes to the formulas in the previous subsection to 
account for special relativity. We replace the suscepti- 
bility in Eq. [8] with the relativistic one of Bergman and 
Eliasson, which is computed using a 3D Jiittner-Synge 
distribution, as opposed to a Maxwellian distribution. 

Given the 3D Jiittner-Synge distribution, 



(Ai = 1.655Ao = 580.88 nm), for which the spatial gain 
rate is = 4 .43 x 10~ 4 ujo/c and ^2/^2 = 0.0465u;o/c. 

We also performed some gain calculations using the 
approximate expression of Bers et. al. for the relativistic 
longitudinal dispersion relation. Their expression shifts 
the gain curve down in wavelength significantly more 
than numerically integrating the formula of Bergman and 
Eliasson. 



fin) 
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the electron susceptibility is given by 



Xe(«2,0 2 ) 
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n d 2 p(fi, K2 /n 2 ) 



K 2 (fi) d/i 2 



(11) 



. (12) 



where 7 is the relativistic factor, /i = m e c 2 /T e , K 2 (fi) 
is the modified Bessel function of the second kind, ^2 = 
u 2 /uo pei anc i ^2 = k 2 c/uj pe . P is given by 



P(fJL,K 2 /Sl 2 ) 



i7r<T . e -M(«2/n2)/V«l/"l-i 



2k 2 /Q; 



(13) 



where a = for > n 2 an d 0" = 1 for f2 2 < ft 2 . 

We can also take into account relativistic effects in the 
electromagnetic dispersion relation, uj 2 = uj 2 e + c 2 fc 2 , by 
using a relativistic version of the plasma frequency: 
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for /i ^> 1. 
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We use this relativistic plasma frequency when we cal- 
culate fco anc i fci- The decrease in the effective plasma 
frequency with temperature is due to relativistic correc- 
tions to the internal energy of the plasma. 31 However, 
this change only has a relatively small impact on the gain 
at the temperature we use in our simulations. It is also 
possible to calculate fco an d fci using the fully relativistic 
transverse dispersion relation from Ref. [28j but we did 
not attempt to do so. 

The overall effect of special relativity on the gain curve 
is shown in Figure [2j Notice that the peak of the relativis- 
tic gain curve lies above the peak of the non-relativistic 
gain curve, and the peak occurs at a shorter wavelength. 
This difference occurs because there are less particles and 
a shallower slope at the plasma wave's phase velocity in 
the Jiittner distribution than in the Maxwellian distri- 
bution, resulting in weaker Landau damping. At the 
peak of the analytic non-relativistic gain curve (Ai = 
1.659A = 582.37 nm), for which k 2 X De = 0.289, the 
non-relativistic damping rate is 3.24 x 10 _3 o;o, while the 
relativistic rate is 2.44 x 10 -3 co>o- The strong damping 
limit still applies at the peak of the relativistic gain curve 



C. PIC Modification 

We can improve the agreement between simulation and 
theoretical results by taking into account a few known 
aspects of finite-difference PIC codes: finite-size par- 
ticles, differencing operators, and field smoothing plus 
compensation.^ The particles have a finite size because 
the charge and current are interpolated to a grid via the 
"shape factor" S(x): 



qS(x) H> qS(x). 



(15) 



To reduce the self-heating and spurious noise from alias- 
ing, we use second-order B-splines. Transforming to 
Fourier space, 



qS(x)^qS(k), 



where 



S(k) 



1 ^sin(fcA/2)Y 
kA/2 ) 



(16) 



(17) 



for ID simulations, with A being the cell width. 

In finite-difference codes, like OSIRIS, differencing op- 
erators modify the dispersion relation by changing the 
relationship between the charge density p, longitudinal 
electric field E 2l and electrostatic potential <j). In Fourier 
space, 



47rp(fc 2 )= k\ 



sin(fc 2 A/2)\ 



a/2 ; 



= K\k 2 )^(k 2 ) 



and 



k 2 A 

= -iK(k 2 )4>(k 2 ). 



(18) 



(19) 



We additionally smooth the fields in our simulations 
to further reduce the effects of aliasing, and we com- 
pensate to reduce numerical modifications to the disper- 
sion relation for small fc. Without the use of splines and 
smoothing, grid heating instabilities occur for A > 3A^ e - 
Finite-difference codes, such as OSIRIS, can use a digi- 
tal filter to compensate for this effect. We perform the 
filtering of some quantity (j) on the grid by replacing 



u with 



1 + 2W 



(20) 
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where j is the grid index and W is a weighting factor. 
Transforming into Fourier space, 



Gain vs. Wavelength 
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1 + 2W 

= SM w (kA)(f> (k)- 



(21) 



We use two filters in the simulations in this paper. The 
first is a two-pass filter, which we use unless stated other- 
wise. The first pass has a stencil of ^ (1,2,1) (W=l/2) and 
the second has a stencil of - (-1,6,-1) (W=-l/6). There- 
fore, 



SM(k) = SM 1/2 (k)SM_ 1/6 (k) 

1 + cos(/cA) 3 - cos(/cA) 



(22) 



The second filter has five passes, and we choose it because 
it causes less deviation from the longitudinal dispersion 
relation without PIC effects than the 2-pass filter. The 
first four passes use a stencil of \ (1,2,1) (W=l/2) and 
the last pass uses a stencil of \ (-5,14,-5) (W=-5/14). We 
perform additional simulations with this filter to observe 
the effect on the gain curve. 



SM(k)= SM* /2 (k)SM_ 
(fcA)\ 



cos 



5/14(*0 
4 



14- lOcos(fcA) 



The particle shape factor, differencing operators, 
field smoothing only affect the plasma frequency, 
simply make the change 



(23) 

and 
We 



k 2 K(k 2 ) 



(L • S(k 2 )) 2 SM(k 2 ) (24) 



everywhere uj pe appears in our formulas to account for 
their effects. 

The effect of the shape factor, differencing operators, 
and both the 2-pass and 5-pass filters on the relativis- 
tic and non-relativistic gain curves is shown in Figure 
[2] Checking the strong damping limit condition for the 
2-pass filter, in the non-relativistic case, the peak gain 
rate drops to g$ r = 2.75 x 10 _4 cjo/c with v 2 jv g2 = 
0.0636cc;o/c, while for the relativistic case, it drops to 
gl = 3.70 x 10- 4 u /c with v 2 jv g2 = 0.0531u;o/c. 



IV. SUBTRACTION TECHNIQUE 

We use a subtraction technique 32 in our data analy- 
sis to clearly see waves with amplitudes below the back- 
ground plasma fluctuation level. The technique requires 
running two simulations, the first with a perturbation 
whose effects we wish to examine, and the second with- 
out the perturbation, but with the same random num- 
ber generator seed. We then subtract the results of the 
second simulation from the results of the first. In our 
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FIG. 2. (Color online) Several calculations of the linear con- 
vective BSRS gain spectra for a box 1790c/cjo (100 /im) long. 
The dash-dotted curves are analytic results, and the (dashed, 
solid) ones take into account PIC effects with the (5,2)-pass 
filter. 



case, the first simulation has both a backward propagat- 
ing light seed pulse and a forward propagating pump, 
while the second simulation has just the pump. 

Figure [3] shows the electrostatic field induced by the 
beating of the pump and the seed pulse. The seed 
pulse in this simulation has a Gaussian-like profile, as 
described in section [V| with Ai = 1.644Ao (577nm) and 
I ls = 5 x 10 -4 io- The amplitude of the plasma wave is so 
small that we cannot distinguish it from the background 
fluctuations without using the subtraction technique. In 
the subtracted result, background fluctuations enter the 
simulation starting at the sides of the box due to the ther- 
malizing boundary conditions. Fluctuations will always 
re-enter the subtracted data as the two simulations be- 
come uncorrelated, but the re-emission of particles with 
random speeds at the boundaries exacerbates this situa- 
tion. 

Because we do not observe meaningful SRS without 
a seed, we use the subtraction technique as a means of 
separating the scattered light from the pump light when 
pump depletion is not significant. When pump depletion 
is small, the subtraction technique for the transverse elec- 
tric field works well for finding the scattered light at all 
positions in the box. However, when pump depletion 
becomes significant, the subtraction technique does not 
produce good results by itself anywhere except at the far 
left side of the box, where pump depletion does not occur. 
Therefore, we can still use the subtraction technique to 
observe the reflectivity at the left side of the box, but we 
need to filter out the pump in Fourier space to observe 
the scattered light anywhere else in the box. 

The number of particles per cell in the simulation af- 
fects the usefulness of the subtraction technique. Fluc- 
tuations take longer to enter the simulation as we in- 
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FIG. 3. (Color online) The longitudinal electric field in a 
simulation with (a) and without (b) a seed pulse. Subtracting 
the latter from the former reveals the plasma wave (c,d). 



crease the number of particles per cell. When we de- 
crease the number of particles per cell, low-intensity seed 
pulses become more difficult to distinguish from the back- 
ground fluctuations when we use the subtraction tech- 
nique, until we cannot distinguish the peak of a pulse 
with Ii s = 5 x 10 -4 io from the fluctuations in simula- 
tions with 512 particles per cell. However, changing the 
number of particles per cell has no significant effect on 
the convective amplification of the seed pulse. 




200 400 600 800 1000 200 400 600 800 1000 
t [1/coJ t [1/co ] 

FIG. 4. The shapes of the seed pulses in the simulations. 
A Gaussian- like pulse (a) with a SOOuJq 1 rise and fall time 
(~ SOOuq 1 FWHM), and a flat-top pulse (b) with a 200^ 1 
Gaussian-like rise and fall, and a steady amplitude for 600^ 1 
in between. 



V. CONVECTIVE GAIN SIMULATIONS & 
COUPLED-MODE RESULTS 

A. OSIRIS Convective Gain Simulations 

We perform OSIRIS PIC simulations to observe seed 
amplification in the linear regime and determine under 
what conditions SRS enters the nonlinear regime. The 
simulations begin at t = with the pump incident from 
the left. After 2,000o;^ 1 , the seed enters from the right, 
by which time the pump has crossed the box. We use 
two different temporal profiles, or "shapes," for the seed 
and vary its wavelength and intensity. The first shape is 
approximately Gaussian, and rises from zero to its peak 
amplitude in r = SOOcj^" 1 , then falls back to zero over an- 
other 500gJo \ for - 500UJQ 1 FWHM. The second shape is 
a fiat-top pulse with a Gaussian-like rise and fall time of 
r = 200co>^" 1 , and a constant peak amplitude for 600^ 1 
in between. The two pulse shapes are plotted in Fig- 
ure [4j We describe the seed pulse using the notation 
I\{z — L, t) = Ii 3 s(i), so that s(t) describes the pulse 
shape and I\ s is the maximum incident intensity. In our 
simulations using seed pulses, I\ s > 5 x 10 -4 io- 

Figure [5] shows the scattered light and the plasma wave 
as a function of position and time in an OSIRIS PIC 
simulation when we use a Gaussian-like seed pulse with 
I l3 = 5 x 10 -4 io and Ai = 1.644A . We also include 
line-outs of the scattered light amplitude vs. position at 
various times to show the evolution of the seed pulse more 
clearly as it crosses the box from right to left. We use 
a Hilbert transform to envelope the results, producing a 
smooth appearance. 

We define the measured gain, g m eas, of a pulse as 



In 



/max(IiO = 0,i)) 

V ha 



(25) 



The plot at the top of Figure [6] shows the gain we measure 
in simulations when we vary the seed wavelength while 
keeping I\ s fixed at 5x 10 _4 io- The linear relativistic gain 
peaks near Ai = 1.644Ao (577nm) in simulations with a 
2-pass filter and near Ai = 1.650Aq (579nm) in simula- 
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FIG. 5. (Color online) Scattered light vs. position and time 

(a) , lineouts of the scattered light vs. position at various times 

(b) , and plasma wave amplitude (c) from a simulation using 
a Gaussian-like seed pulse with Ii s = 5 x 10 _4 /o and Ai = 
1.644A . 



tions with a 5-pass filter. Our simulation results agree 
with these predictions. We also plot a gain curve from 
simulations using a non-relativistic particle pusher and 
a Maxwellian velocity distribution. This non-relativistic 
curve lies below the relativistic one at most points, as we 
expect, and the location of its peak agrees with the gain 
curve from non-relativistic theory with a 2-pass filter. 
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FIG. 6. (Color online) Measured gain in simulations with box 
length 1790c/cjo (100 fim) vs. seed wavelength (a) and seed 
intensity ii s (b). Plotted are the gains we measure using a 
Gaussian-like pulse (green with square markers) and a flat- 
top pulse (blue with circle markers), both with Ai = 1.644Ao- 
(a) includes the gain from simulations with a Gaussian-like 
pulse and the 5-pass filter (black with 'x' markers), plus those 
without relativistic effects (brown with '+' markers). Several 
theoretical gain curves taking into account PIC effects are 
included: 2-pass non-relativistic is dash-dotted red, 2-pass 
relativistic is solid magenta, and 5-pass relativistic is dashed 
purple. 



In the plot on the bottom of Figure [6j we see how the 
measured gain changes as we vary 7 ls while keeping the 
seed wavelength fixed at Ai = 1.644Ao- The measured 
gain of the Gaussian-like pulse remains relatively con- 
stant as we increase the initial seed intensity, until the 
seed intensity reaches several times the pump intensity. 
This behavior indicates that we are in the linear regime 
for even large seed amplitudes. The deviation from linear 
theory at the highest amplitudes is due to pump deple- 
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tion. However, the measured gain of the flat-top pulse 
increases with seed intensity before falling off. This in- 
flationary gain is caused by the higher amplitude seed 
pulse generating a larger amplitude plasma wave and, in 
the presence of the seed pulse, particles executing several 
bounces, thus decreasing the Landau damping rate. We 
explain this effect in more detail in the following sections. 
We note that if we had decreased (increased) the seed 
pulse length, the deviation from linear behavior would 
occur at higher (lower) seed intensity. 



Reflected Light, Flat-Top Seed 




B. The Coupled- Mode Equations 

The measured gain of a seed pulse can differ from the 
steady-state linear result due to several linear and non- 
linear effects. Linear effects include pulse shape, with 
each frequency in a spectrum of incident frequencies be- 
ing amplified at a different rate. Nonlinear effects include 
pump depletion and nonlinear (kinetic) and nonlocal re- 
ductions to the real part of the frequency and the damp- 
ing rate of the plasma wave. In this section, we investi- 
gate the effect of the pulse shape on th e measured gain 
using the coupled- mode equations! 33 * 34 ! A comparison of 
the coupled-mode and OSIRIS results, similar to that 
performed by Wang et. al. in Ref. [35j isolates linear from 
truly nonlinear, kinetic physics, and provides confidence 
in the PIC method. 

In the coupled- mode equations, we let (c^, hi) of the 
carrier waves be real and work with complex envelopes 
di(x,t). We assume that the envelopes vary slowly with 
respect to the carriers, such that |Va^| <C |fc^a^| and 
\ddi/dt\ <C |c<^a;|. The complex envelopes for the action 
amplitudes dj relate to the physical quantities by 



Aa = -i 



2tt 



1/2 



i j — o i — i cij ex.p[i(kj • x — ujji)]k + cc, j = 0, 1 



(26) 



for light waves, with Aj the vector potential, and by 

n\ = -— a 2 expk(/c 2 • x - cj 2 £) + cc, (27) 

2 \m e U2 J 

for the plasma wave. n# is the spatially varying back- 
ground electron density, n\ is the perturbation on top of 
it, and cc denotes complex conjugate. 
The couple mode equations are 

d \ 

-Q t + v g o • V + i/ + i8q J a>o= Kaia 2 , (28) 

d_ 
dt 
d_ 

at 



v g i ' V + v\ + iSij di= —Kd^d\, (29) 

Vg2 • V + V2 + i$2 I Q>2= ~KdQd\, (30) 



where the coupling constant is 



K 



FIG. 7. (Color online) The BSRS reflected light using a flat- 
top seed with Ai = 1.644Ao (a) and 1.650Ao (b). The simu- 
lation results are plotted with solid magenta lines while the 
coupled-mode results are plotted with dash-dotted black lines. 
The horizontal dashed red line indicates the maximum seed 
amplitude with no gain. 



Vi is the damping rate of mode z, and Si = (ujp e — 
uJf + c 2 kf)/(2uji) is the detuning frequency for the light 
waves, which reflects the departure of mode i from be- 
ing a natural mode of the plasma. We assume no in- 
verse bremsstrahlung, which is the case in the PIC sim- 
ulations, so v$ = v\ = 0. The Landau damping rate, 
^2 = Si/(de r /du2), where e r = Re[e], S{ = Im[e], and e 
is the kinetic dielectric function. We set the light- wave 
detuning, So = Si = 0. For the electrostatic mode, we 
find S2 using the kinetic equation, 5 2 = —e r /(de r /duj2). 



^ujqujiuj2 \/8riBm e ' 



(31) 



The equations in Sections III B and III C allow us to take 
into account special relativity and PIC effects when we 
calculate the coefficients. 

We can directly compare the results from OSIRIS and 
the coupled- mode solver by examining the reflected light. 
Figure [7] shows the reflected light from runs with 7 ls = 
5 x 10 io using flat-top pulses with Ai = 1.644Ao and 
1.650Ao- In Figure [8j we compare the measured gain from 
PIC simulations and the coupled-mode solver for various 
wavelengths. The simulation and coupled-mode results 
are in excellent agreement for both the Gaussian-like and 
flat-top pulse runs. 

The agreement between simulation and coupled-mode 
results is not as good when we examine the longitudinal 
field. Figure [9] shows the amplitude of the longitudinal 
field at t = 4,000o;^ 1 for the same runs as shown in 
Figure [7) The disagreement occurs soon after the wave 
begins growing and is visible at about x = 300c/ 'ujq. We 
are not yet sure of the reason for this disagreement. The 
disagreement becomes worse as background fluctuations 
begin to enter the PIC results after x = 400c/cjq. 



VI. INFLATION AFTER SEED PASSAGE 

For short seed pulses, we find significant reflectivity 
well after the seed propagates out of the simulation box. 
The pump now directly interacts with the Langmuir wave 
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FIG. 8. (Color online) The measured gain in a box of length 
1790c/cjo (100 jLtm) predicted by the coupled- mode solver as 
we vary the wavelength using a Gaussian-like pulse (dashed 
brown with markers) and a flat-top pulse (dashed ma- 
genta with circle markers). For comparison, we include the 
PIC simulation result for the Gaussian-like pulse (solid blue 
with c x' markers) and the flat-top pulse (solid red with square 
markers) along with the linear relativistic gain curve taking 
into account PIC effects (dash-dotted green). 
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FIG. 9. (Color online) The longitudinal field vs. space at 
t — 4, OOOcj^ -1 when using a flat-top seed with Ai = 1.644Ao 
(a) and 1.650Ao (b); the same runs as used in Fig. [7] The 
simulation results are plotted with solid magenta lines while 
the coupled-mode results are plotted with dash-dotted black 
lines. 



that is still present in the plasma after the seed leaves 
the box. Without trapped particles, the Langmuir wave 
is described by its linear dispersion relation, BSRS re- 
mains in the strongly-damped limit, and no observable 
growth of BSRS occurs for the pump intensity and the 
plasma length of interest. However, a small amplitude 
wave can evolve into a nonlinear weakly dampe d wav e af- 
ter the trapped particles execute a few bounces! 14 * 15 * The 
period for a bounce, or bounce time, is tb = 2tt/ujb, 
where ujb = y eE2k/m e is the bounce frequency for 
deeply trapped electrons, E<i is the electric field ampli- 
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FIG. 10. (Color online) The scattered light (a) and the plasma 
wave (b) seen when we extend the duration of the simulation 
using a Gaussian-like seed pulse with l\ s = 8 X 10 _3 /o and 
Ai = 1.644Aq. The seed exits the box around t — 5, OOOcj^ -1 . 



tude, and k is the wavenumber of the wave. As particles 
are trapped, the damping rate decreases below its lin- 
ear value to a residual level which depends on details of 
the problem!^ The ponderomotive beating of the pump 
and the scattered light will drive the wave to increasing 
amplitudes. Such a situation will lead to noticeable re- 
flectivity later in the simulation. The seed pulse must be 
the cause of any such reflectivity because, with the pump 
amplitude we use in our simulations, BSRS is negligible 
without a seed. BSRS that occurs after the seed has 
left the box is both useful for isolating the process of ki- 
netic inflation^ and is potentially relevant to situations 
where BSRS in one region of space or time seeds BSRS 
in another one, creating a plasma wave that triggers an 
inflationary burst of BSRS. 

Figure IlO] shows the scattered light and plasma wave 
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FIG. 11. (Color online) Evolution of the scattered light (a) 
and the longitudinal field (the PIC simulation field is filtered) 
(b) with time at x = 550c/cjo for the simulation of Figure [To] 
The simulation results are plotted with solid blue lines and 
the coupled-mode results are plotted with dashed black lines. 



for times after the seed leaves the box in a simulation 
using a Gaussian-like seed pulse with I\ s = 8 x 10 -3 io 
and Ai = 1.644Ao- We observe high reflectivity after 
t = IC^OOOcj^" 1 , along with a corresponding growth in 



the plasma wave. Figure 11 shows lineouts of the scat- 
tered light and longitudinal field at x = 550c/co>o along 
with a comparison with the coupled- mode result. For 
the longitudinal field from the PIC simulation, we filter 
out all modes except IAujq/c < k < l.bujo/c. Notice the 
dip in the plasma wave amplitude around t = 4, 500^ x , 
corresponding to the drop in the seed's amplitude, before 
the plasma wave begins to grow again. 

We verify that kinetic inflation is occurring by tracking 
particles traveling near the plasma wave phase velocity 
and plotting their orbits in the frame of the wave, as done 
in Figure [12] from t = 5, OOOcJq" 1 to t = 7, OOOcJq" 1 . The el- 
liptic trajectories are clear indicators of particle trapping. 
During this time, the distribution function begins to flat- 
ten around the Langmuir wave phase velocity, as seen in 



Figure[13l which is another indication of particle trapping 
and a clear indication of the reduction of Landau damp- 
ing. The tail is flattened to much higher velocities during 
the larger burst of SRS that grows after the seed leaves; 
for example, as shown at t = 17, OOOcj^" 1 . This larger tail 
in our simulations potentially contributes to the energetic 
electrons (commonly referred to as "hot electrons") seen 
in recent experiments.^ The production o f hot electron 
tails by SRS is an active area of researchf 38 * 39 1 but the 
small flattening early in time is sufficient to affect the 
growth and further onset of SRS we are studying here. 
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FIG. 12. (Color online) Orbits of particles trapped in the 
plasma wave in the simulation of Figure |1Q| The orbits are 
plotted for 2,000cj^" 1 starting at to = S^OOOuq 1 in the wave 
frame, where v p h = 0.271c is the plasma wave phase speed and 
the orbits are centered around the relativist ic p p h = 0.282m e c. 



Measured and Juttner Distributions 
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FIG. 13. (Color onli ne) The distribution functions from the 
simulation of Figure 
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at t = GjOOOcJq -1 (solid blue) and 
t — 17, OOOcJq -1 (dashed black) along with the Juttner dis- 
tribution for a 2.5keV electron plasma (dash-dotted red). 
The measured distributions have a plateaus beginning around 
7/3 = 0.26, indicating particle trapping. 



Since particle trapping causes the inflationary scatter- 
ing at late times, the inflationary bursts will occur earlier 
if the seed drives a larger Langmuir wave with a shorter 
r#, so that the trapped particles accumulate bounces 
faster. We can increase the Langmuir wave amplitude 
by increasing the intensity, the duration, or choosing a 
seed wavelength that produces a higher gain. 

We first vary the initial amplitude of the seed. Figure 
[14| shows that the first burst of reflected light has a maxi- 
mum around t = 17, OOOa;^ 1 when we use a Gaussian-like 
seed pulse with I ls = 4x 10 _3 / and Ai = 1.644A . As we 
increase ii s , the burst moves earlier, but the difference 
in amplitude between the burst and final seed amplitude 
also decreases. As a clear demonstration of the effect 
of the higher intensity seeds, in Table [I] we examine the 
plasma wave amplitudes and particle bounce times for 

We examine a selection of 



the simulations in Figure 14 



particles near the plasma wave phase velocity between 
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I la Avg. EPW Amp. 

M [E ] 



Calc. 



TB 



Meas. tb 



4 x 10" 
8 x 10" 
0.016 
0.032 



6 0.020 
3 0.033 
0.049 
0.071 



610 
480 
390 
320 



675 
550 
415 
350 



The Wigner transform results for a run with a flat-top 
seed with I ls = 8 x 10 -3 2o and Ai = 1.638Aq are shown 



TABLE I. The bounce times of deeply trapped electrons mea- 
sured in the simulations of Figure |14[ along with the bounce 
times calculated using the average plasma wave field ampli- 
tude along the particle's trajectory. 



Reflected Light, Extended Duration Simulations 




FIG. 14. (Color online) The reflected light in extended 
duration simulations using a Gaussian-like seed pulse with 
Ai = 1.644Ao and various maximum initial intensities. Shown 
are I ls = 4 x 10~ 3 / (solid blue), 8 x 10~ 3 / (dotted green), 
0.016 Jo (dash-dotted red), and 0.032/ (dashed black). 



x = 350c/co>o and x — 450c/co>o at t = b^OOOuj^ 1 . We use 
the most deeply trapped particle to determine the time 
it takes to complete one bounce. To calculate the theo- 
retical bounce time, we filter out all plasma wave modes 
except IAujq/c < k < 1.5ojq/c and average the field am- 
plitude over the bounce time along the particle's track. 
We then substitute the measured values of the amplitude 
and wavenumber (k = IAAujq/c) into the formula for the 
bounce time. The results are in Table U The measured 
bounce times are slightly (~ 10%) longer than the simple 
expression. We believe this discrepancy is due to the fact 
that the amplitude of the wave is changing with time and 
because the calculation is for a parabolic potential well, 
while the particles are actually trapped in a sinusoidal 
well. 

Besides lowering the kinetic damping rate, another 
well-known effect of particle trapping is the nonlinear fre- 
quency down-shift of the plasma waveP^ As the plasma 
wave grows, it will shift downward in frequency because 
it will trap more particles. According to the frequency 
matching condition in Eq. [Ta| the down-shift in the fre- 
quency of the plasma wave should be accompanied by 
an up-shift in the frequency of the scattered light. We 
examine this down-shift using a Wigner transform with 
a Choi filter The Wigner transform takes a function 
of time and computes its representation as a function of 
both frequency and time. It maps f(t) — >• f(uj,t). 



in Figure 15 The results for the same run, except using a 
seed with Ai = 1.627Aq, are in Figure 16 In both cases, 



the seed appears in the Wigner transform scattered light 
plots around t = 4500o;^ 1 , when the seed reaches the left 
side of the box. The scattered light frequency shifts up 
while the plasma wave frequency shifts down, as expected 
due to trapping. 

The inflationary bursts begin growing near the fre- 
quency with the highest gain regardless of the seed's cen- 
tral frequency, consistent with a harmonic oscillator that 
is driven off-resonance. The initial growth is near the cen- 
tral frequency of the seed in the case with Ai = 1.638Ao, 
and in the seed's lower-frequency tail in the case with 
Ai = 1.627Ao- This observation suggests that the pon- 
deromotive beating of the seed and the pump simply dis- 
turbs the plasma and provides an initial level for growth, 
but the subsequent growth occurs at the most unstable 
mode. 

There are two inflationary bursts in the scattered light 
plot in Figure [15} According to Ref. [T71 the separation 
of the bursts in time is 2tt/Aujnl, where Auonl is the 
nonlinear shift of the scattered light and plasma wave 
from their non-inflationary resonant frequencies, where 
peak linear gain occurs. In our case, the bursts occur at 
different frequencies. The first is roughly on-resonance 
while the second has Auj^l ~ 0.002co>o. Averaging the 
two, we have Auj^I = O.OOIcjo, or a separation of about 
6,000o;^ 1 , which is in good agreement with the actual 
separation of the two bursts. Figures [15] and [16] show 
the seed wavelength also affects the time required for 
inflation to set in. We explore this effect in Sections |VII| 
and [Villi 

Increasing the width of the seed drives the plasma wave 
longer, increasing its amplitude and causing inflationary 
bursts to occur earlier. In Figure [l7| we vary the duration 
of a flat-top pulse while keeping the rise and fall time 
at 200u;^ 1 . The burst of high reflectivity moves earlier 
and, correspondingly, the plasma wave reaches a high 
amplitude quicker as we increase the pulse duration from 
^OOOco^ 1 to 3,000a;^ . However, this effect only occurs 
if we drive near resonance. If we drive off resonance, 
then the the large bursts don't occur, and the reflectivity 
oscillates with a period of 2tt/Aujnr, where Aunr is the 
difference between the seed frequency and the resonant 
frequency, at which inflationary SRS grows. 



VII. INFLATION OF THE SEED 

When the seed pulse width becomes comparable to the 
bounce time, then the seed itself can undergo inflation, 
illustrating the difference between the linear and infla- 
tionary regimes. We perform simulations using a flat-top 
pulse with a duration of 6, OOOuJq 1 with / ls = 8 x 10 -3 io 
at various wavelengths. The reflectivity plot in Fig- 



ure 18 demonstrates the variation of the scattered light 
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FIG. 15. (Color online) The scattered light (a) and plasma 
wave (b) in a simulation using a flat-top seed pulse with Ii s = 
8 x 10" 3 / and Ai = 1.638A . Below them are the Wigner 
transforms of the reflected light at x = (c), the plasma wave 
at x = 500c/cjo (d), and the plasma wave at x = 900c/o;o (e). 



with time in the simulations as compared to the steady- 
state values from linear theory. The reflectivities seen in 
the simulations can significantly exceed the linear the- 
ory values, particularly for the seeds with Ai = 1.638Ao 
and 1.644Ao- The runs using seeds with Ai = 1.632 An 
and 1.650Ao reach levels above the linear theory values, 
but dip below linear values several times due to driving 
off resonance. In these four cases, inflationary scatter- 
ing continues after the seed ends around uo^t — 10, 200. 
When we use seeds with Ai = 1.627An and 1.658 An, the 
reflectivity does not reach significantly above the linear 
value and drops down when the seed ends, because infla- 
tionary scattering does not occur easily when the seed is 
far from resonance. 

We observe oscillations in the reflectivity in all the 
cases, except when we use the seeds with Ai = 1.638An 
and 1.644An- These oscillations are due to the scattered 



20000 



15000 



100001 



- 50001 




1UUU 
X [C/0) Q ] 





0.618" 


-4 .. 


^0.616: 


ude 


0.614- 


< 

-6 


| 0.612 1 
1" 0.610: 
0.608^ 





Central Seed Frequency 


L 





1500 




10" 3 




10" 4 -g 




io" 5 i_ 








10" 6 < 




10" 7 



? 0.390 
. 0.388 
j 0.386 
[ 0.384 
0.382 



time [l/co ] 

Wigner Transform of the Longitudinal Field 



(d) 



X=900C/G) n 







5.0 10 3 



0020 

.0015 % 

.0010 S. 
E 

.0005 < 



1.0 10^ 
time [l/co ] 



1.5 10^ 



0.0000 



FIG. 16. (Color online) The scattered light (a) and plasma 
wave (b) in a simulation using a flat-top seed pulse with I\ s — 
8 x 10~ 3 / and Ai = 1.627A . Below them are the Wigner 
transforms of the reflected light at x — (c) and the plasma 
wave at x — 900c/cjo (d). 



light seed driving SRS off resonance.^ The ponderomo- 
tive beating of the pump and scattered light drives the 
plasma wave, so if the phase of the beat drive leads the 
plasma wave by 7r/2, the plasma wave no longer grows. 



The plasma wave density ri\ 



drive F p oc 



dx 

-^EqE\. The product, 



£?2, while the beat 



dx 



dx 



E E 1 



(32) 



indicates the phase difference between the beat drive and 
the plasma waveP^ We also note that the change in the 
pump energy density with time is given by 



dW 



87rm e co>i/c2 ^ P 



(33) 



Therefore, if R p is positive, the waves are in phase, the 
beat drive is resonantly driving the plasma wave, and en- 
ergy is transferred from the pump to the seed and plasma 
wave. The inverse applies if R p is negative. 

Figure [19] shows the scattered light and the resonance 
plot for the seed with Ai = 1.644An- We smooth the re- 
sult from the resonance diagnostic in x using a 6-point 
moving average. Notice the dip (valley) in the scattered 
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FIG. 17. (Color online) Results from simulations with flat- 
top seed pulses with I\ s — 8 x 10 _3 /o and Ai = 1.644Ao- The 
scattered light (left) and plasma wave (right) for simulations 
using a pulse duriation of (1000, 3000)0;,^ 1 are in panels (a-b, 
c-d), respectively. 



Reflected Light, Extended Seeds 
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FIG. 18. (Color online) The reflected light seen in several 
simulations using a flat-top seed of duration 6, OOOuJq 1 with 
hs = 8 x 10 _3 /o for various wavelengths. The two top curves 
are for seeds with Ai = 1.638Ao (solid red) and 1.644Ao 
(dashed cyan). The four lower curves are for seeds with 
Ai = 1.627A (dotted blue), 1.632A (dashed green), 1.650A 
(dash-dotted purple), and 1.658Ao (solid yellow). For compar- 
ison, we mark the steady-state linear relativistic PIC values 
using horizontal dashes on the left side of the plot. The red 
vertical dash on the lower right side of the plot indicates ap- 
proximately when the seeds end. 



light amplitude around x = 1250c/u;o, t = 7, OOOu;^ 1 , 
and the corresponding negative area on the resonance 
plot around x = 1500c/co>o. When this drop occurs, the 
beat drive and the plasma wave are out of phase, so en- 
ergy flows from the scattered light wave to the pump. 
This shift away from resonance is due to the nonlinear 
frequency shift of the plasma wave, explained in Section 
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FIG. 19. The scattered light (a), and the resonance product 
R p from Eq. 32 (b) for the flat-top seed of duration 6, OOOuJq 1 
with Ai = 1.644A . 



VI The non-resonant drive is also responsible for the os- 



cillations we see in the reflectivity plot of Figure 18 



Figure [20] shows the same plots as Figure [19j except 
using a seed with Ai = 1.638Ao- This seed continuously 
drives the scattered light to its maximum amplitude be- 
fore dropping off, whereas the one with Ai = 1.644Ao 
drives a lower-amplitude burst of scattered light before 
it produces a second burst at much higher amplitude. 



Notice that the resonance plot in Figure 20 shows that 
the beat wave and plasma wave are in resonance until 
the peak of the scattered light burst. This resonant drive 
leads to a burst of scattered light at higher amplitude 
than in the Ai = 1.644Aq case. 
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FIG. 20. The scattered light (a), and the resonance product 
R p from Eq. 32 (b) for the flat-top seed of duration 6, OOOuJq 1 
with Ai = 1.638Aq. 



VIII. INFLATION OF CONTINUOUS SEEDS 

In this section, we extend the flat-top seed pulse so that 
it remains on through the end of the simulation. Based 
on what we have observed in the previous sections, we 
expect the inflationary behavior we see in these simula- 
tions to depend on the seed intensity and wavelength. 
We know that there must be an intensity threshold be- 
low which the seed does not drive inflationary behavior 
in the simulations (very long bounce times), because we 
see negligible scattering without a seed. However, as we 
saw in Section |YlI[ a seed that is intense enough to drive 
inflation on resonance may not be intense enough to drive 
it off resonance. 

Furthermore, the use of continuous seeds is related to 
past work by Winjum et. al. on scattering off of plasma 
wave packets. 17 Light scattering off a wave packet that 




'is V 

FIG. 21. (Color online) The time-average measured gain seen 
in simulations as we vary the continuous seed intensity using 
seeds with Ai = 1.644Ao (dashed blue) and 1.658Ao (solid 
red). For comparison, we use horizontal dashes on the left 
side of the plot to mark the steady-state gain from linear 
relativistic PIC theory. 



has undergone a nonlinear frequency shift acts as a seed 
for BSRS in the unperturbed background plasma. How- 
ever, the beat drive frequency is not a natural mode of 
the background plasma. This non-resonant beat drive 
leads to oscillations in the reflected light with a period 
2tt/Aujnl, where Aujnl is the nonlinear frequency shift 
of the plasma wave, much like the amplitude of a simple 
harmonic oscillator varies when driven off resonance. In 
this section, we examine the effect of resonant and non- 
resonant drive in more detail by using continuous seeds 
at different wavelengths and intensities, and examining 
the results using plots of the reflected light. 



Figure 21 shows the time- average measured gain for 
various seed intensities with Ai = 1.644Ao and 1.658Ar> 
When we use a seed with Ai = 1.644Ao, an intensity 
of 1.25 x 10 -4 io (smallest value shown) is enough to 
cause inflation. However, when we use a seed with 
Ai = 1.658Ao, we do not see inflation until the seed in- 
tensity reaches 8 x 10 -3 io- The measured gain at both 
wavelengths decreases with seed intensity due to pump 
depletion once inflation sets in. 

In order to understand the results with a 1.658Ao seed, 
we first discuss trapping effects on the SRS gain spec- 
trum via Eq. 10 Trapping nonlinearity leads to a re- 



duction in Xi an d thus Landau damping, as well as a 
down-shift in the natural plasma wave frequency. This 
latter effect decreases the resonant scattered wavelength, 
where 1 + Xr = 0. In electrostatic simulations with an 
external driver, Fahlen showed behavior consistent with 
this picture.^ For feA^e ~ 0.3 and a driver frequency 
above the natural frequency, the plasma response is even 
smaller than the linear response. However, when the 
driver is below the natural frequency, a larger response 
is obtained. 
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Trapping, therefore, makes the gain spectrum nar- 
rower and peaked at a smaller wavelength. The gain 
increases for wavelengths near resonance. However, for 
wavelengths far from resonance, the gain does not in- 
crease, and, in fact, vanishes as \% ~^ 0- We can see this 



effect clearly by examining Eq. 10 Since the seed is far 
from resonance, 1 + Xr ^> Xu so 9o oc Xi/0- + Xr) 2 - This 



vanishing gain is clear in Figure 22 1, where the gain ap- 
proaches zero in steady state for the continuous 1.658Ao 
seed with h s /Io = 4 x 10 -3 . Thus, inflation is not pos- 
sible at non-resonant wavelengths. For hot, low-density 
plasmas, there is no resonant wavelength, and all phase- 
matched plasma waves satisfy the loss of resonance con- 
dition feADe > 0.53.^ Inflation cannot occur in such a 
plasma at any scattered wavelength. 



The red curve in Figure 21 is for a seed wavelength that 
is non-resonant and larger than the linear resonance. The 
reduction of Xi due to trapping cannot lead to inflation 
at this wavelength, and the nonlinear frequency shift will 
move the resonance farther away. Both effects conspire to 
reduce the SRS gain below its linear value, which is what 
we observe for the lowest seed intensities. The num ber of 
bounce orbits completed by resonant electrons j 13 * 43 * based 
on the plasma wave amplitude computed from linear the- 
ory, is > 4. It is thus consistent for trapping nonlinear ity 
to occur and reduce the SRS gain. The increase in gain 
for Iis/Iq = 8 x 10 -3 , as shown in Figure 22 3, first de- 



velops at linear resonance, 1.644Ao, not at the seed value 
of 1.658Ao, then shifts down in wavelength with time to 
finish around 1.638Aq. This progression is similar to the 



one in Figure 16 We performed a similar run using a 
fiat-top seed of duration I^OOOuJq 1 with a central wave- 
length of 1.658Ao and hs/Io = 4 x 10 -3 and observed 
inflation similar to that in Figure [22] p, but without the 
oscillations, which raises the possibility that continuous 
seeds can suppress inflation. 

There are oscillations in Figure [22) 3, which occur be- 
cause the seed is driving SRS off resonance. As we dis- 
cussed earlier, this non-resonant drive leads to oscilla- 
tions in the reflected light with a period of 2tt/Aujnr, 
where Aujnr is the difference between the seed frequency 
and the resonant frequency, at which inflationary SRS 
grows. Equivalently, the Aunr is the difference be- 
tween the seeded beat drive frequency and the frequency 
of the plasma wave packet. When we use a seed with 
hs/^o = 1.024, the oscillations are more prominent and 



faster, and we see no amplification, as seen in Figure 22 2. 
The increased oscillation frequency is due to the higher 
amplitude plasma wave undergoing a greater frequency 
shift. We also see a beat wave pattern covering many 
oscillations, which is caused by the nonlinear frequency 
shift of the plasma wave packet, as described earlier in 
this section. 



Figure 23 shows the reflected light in simulations using 
a (low-, moderate-, high-) intensity seed with hs/Io = 
(1.25 x 10- 4 ,8 x 10- 3 , 1.024), but Ai = 1.644A . When 
we use the low-intensity seed, we see the reflected light 
increase monotonically until it saturates near the end of 
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FIG. 22. (Color online) The reflected light in simulations 
using a continuous seed with Ai = 1.658Ao (solid red in Fig. 
2lJ. I ls /Io = 4 x 1(T 3 (a), 8 x 1(T 3 (b), and 1.024J (c). 
The horizontal dash-dotted line indicates the amplitude of 
the seed (unamplifled) reflected light. 
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FIG. 23. (Color online) The reflected light in simulations 
usin g a continuous seed with Ai = 1.644Ao (dashed blue in 
Fig. |2lJ. I ls /Io = 1.25 x 10" 4 (a), 8 x 10" 3 (b), and 1.024 
(c). The horizontal dash-dotted line indicates the amplitude 
of the seed (unamplifled) reflected light. 



the simulation. Unlike when we used Ai = 1.658Ao, this 
seed is near resonance, so 1 + Xr << Xi-> an d go cx 1/xi- 
Therefore, particle trapping increases the gain, as ex- 
plained earlier in the paper. Oscillations begin to appear 
again in the simulation with the moderate-intensity seed, 
which indicates that the seed is slightly off resonance, 
as we expect to occur as the plasma wave undergoes a 
nonlinear frequency shift. These oscillations once again 
become more prominent when we use the high-intensity 
seed, and we see a beat wave pattern appear again. 
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IX. CONCLUSIONS AND FUTURE RESEARCH 

Using ID OSIRIS PIC simulations, we have studied 
BSRS of a well-defined seed pulse with variable shape, 
intensity, and wavelength. We found that backward Ra- 
man amplification of a seed scattered light pulse can re- 
main in the strongly damped convective regime. That 
is, for a sufficiently weak seed, kinetic inflation does not 
occur. Peak seed amplification occurs near the peak of 
the linear gain spectrum when we take into account spe- 
cial relativity and PIC effects such as finite particle size, 
finite- difference operators, and field smoothing. Detailed 
comparisons with linear coupled-mode predictions for en- 
velope dynamics show excellent agreement. If the seed 
pulse is intense enough, the driven plasma wave traps 
particles, thereby lowering the Landau damping after a 
time on the order of a bounce period. The plasma wave 
continues to scatter light and grows after the seed leaves 
the system, leading to kinetic inflation later in the simula- 
tion. When we extend the seed in time, kinetic inflation 
occurs while the seed is still present, and we measure 
dynamic seed amplification, which significantly exceeds 
the linear gain rate at times, and can also turn negative 
at other times (i.e., it transfers energy into the pump). 
When we use a continuous seed, we find that the onset of 
inflation depends strongly on the seed wavelength. When 
the continuous seed wavelength is near the peak of the 
gain curve, we see kinetic inflation occur with very low 
seed amplitudes, while higher seed amplitudes are neces- 
sary to drive kinetic inflation using a non-resonant seed. 

Connecting kinetic inflation with experiments involv- 
ing lasers with many hot spots or speckles is still at an 
early stage. Past research has demonstrated that hot 
electrons, beam acoustic modes, and side-scattered light 
can couple hot-spots transversely in 2D simulations, lead- 
ing to higher BSRS from each hot spot than one would 
see without coupling! 19 * 44 * Our results suggest that we 
can expect to see inflationary scattering from the lower- 
intensity parts of the beam interacting with flattened 
(non-Maxwellian) distributions or plasma waves gener- 
ated in more intense parts, or increased scattered light 
from a few speckles triggering a "chain reaction" of down- 
stream inflation. If this chain reaction occurs, most of 
the Raman in underdense laser-produced plasmas, such 
as ICF targets, will be inflationary. PIC codes can be 
used to model hundreds of speckles, but to simulate large 
volumes across a hohlraum effectively requires the use of 
envelope codes (such as pF3D), for which reduced mod- 
els of kinetic nonlinearity are being pursued by several 
groupsPEZl 

An important factor we have not explored is de- 
trapping mechanisms. All of our simulations in this pa- 
per are ID and "collisionless," which means that plasma 
waves easily trap electrons, and trapped electrons cannot 
leave sideways or be scattered out of the plasma wave. 
Since particle trapping is necessary for inflationary scat- 
tering, any de-trapping mechanisms make it less likely 
to occur and impose a threshold amplitude for inflation 



(the threshold in the present work is set by the finite 
amplitude and duration of the seed). In particular, colli- 
sions can kick electrons out of a plasma wave's potential 
well, and electrons can traverse the plasma wave in less 
than one bounce period in higher dimensions! 43 * 48 * Future 
research in this area should explore the effect of these 
de-trapping processes on inflation, and validate reduced 
descriptions. 
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